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1. Introduction. The use of logarithms is so general that it seems 
desirable to give here a short account of the history of their invention and of 
the best method for their computation. To understand these explanations is 
necessary a comprehension of the following definitions and elementary proper- 
ties of logarithms : 

If a" = JV, where a and iV" are positive numbers, x is called the "loga- 
rithm of iV with the base a." We write the relation in the form 

X = log„ N'. 

Two bases are in common use, the base ten used in numerical work, and 
the "natural" base, usually denoted by e, used in the calculus. The value of e 
is the sum of tlie infinite series : 

1 1 1 1 

^ = l + l + 2! + 3!+--- 

= 2.71828 

(1) 
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correct to five decimals. If M is any positive number, 

log„ iV+ logo -M'= logo J^M- 
lip is any real number, 

logoiV>=j)log„J\r. 
From these equations it follows that 



, T.^ logs ^ 1 

log„ N = -p2£ — = 

^" log^a logj^o' 



where b is any positive number. If then the logarithms with one base b, of 
all positive numbers are known, the logarithms of all positive numbers with 
any base a, may be found by division. 

We shall hereafter denote logjo iV"by log iV, and loge N by IN. With 
this notation we may write 

It is proved in the calculus that 

(1) Z(l + x)=x_^+|-^+... if|a;|<l, 

where the symbol | x \ means the absolute or numerical value of x. From ( 1 ) 
it follows that, if | a; [ be less than one, 

(!') log„(l +a;) =log„e |x-|- + --^+ ...j. 

The factor, logo c, is called the modulus of the system of logarithms with the 
base a. We shall denote by ilfthe modulus of the system of base ten. The 
value of -3f, correct to five decimals, is 

lf = log 6 = 0.43429. 

The series ( 1 ) is not convenient for the computation of logarithms, first 
because it does not converge if | a; ( is greater than one, secondly because if 
I a3 1 is only a little less than one a very large number of terms must be taken 
to obtain an approximately correct value for Z(l + a;). But the series is of 
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great importance in any discussion of the properties of logarithms, and more- 
over serves as a basis for the deduction of other series of great use in com- 
putation. One of the most useful of such series is the following : 

i}/ = ^i(y+i) + i i{y - 1) 
(2) 

1 1/ 1 \» 1 / 1 \» 

"^ 2y« - 1 ■*" 3 by** - i) "*" 5 V2y*- l) + ■ ■ ■ • 

This series converges for all values of y greater than one, and converges very 
rapidly for large values. 

It may be deduced from (1) as follows : Replacing x by — a; in (1), 

0^2 S*^ ^* 

Z(l-a;) ^-a:-- -- _- |a;| < 1- 

Subtracting this series from (1) and dividing by two, 

1 , 1 -I- a; 7? X* 

Let X be written 

1 

x = 



Then in order that | a; | be less than one it is both necessary and sufficient that 
y* be greater than one. We find 

1 1 / 1 y 1 / 1 V 

from which series (2) is at once obtained. 

For logarithms with the base ten, we have 

log y = J log (y + 1) + J log (y - 1) 
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If y is a prime integer greater than two, the nunabers y + 1 and y — 1 are 
both even, and their logarithms may be expressed in terms of the logarithms 
of two and of integers less than y. 

2- Historical. A detailed and interesting account of the invention of 
logarithms and of the methods of the first computers is given in the article on 
"logarithms" by J. W. L. Glaisher in the Encyclopedia Brittanica. From 
this I have taken most of the historical matter here given, and to this I refer 
the reader who desires a fuller account than is contained in this article. 

Logarithms were invented, it is generally admitted, about 1614 by an 
Englishman, John Napier, Baron of Merchiston. His first published work 
was seen by Henry Briggs, at that time Professor of Geometry at Gresbam 
College, London, later a professor at Oxford. Briggs visited Napier and 
worked with him towards the perfection of the theory, and subsequently de- 
voted much time to the calculation of logarithms. From these two men the 
system of logarithms with the base ten has derived the name of Naperian or 
Briggsian logarithms. The word "logarithm," it may here be explained, 
comes from the Greek words, \6ya>v api6fio<i, meaning the number of ratios, 
for logarithms were first regarded as a number of ratios. Thus if ten be re- 
garded as the product of 10000 equal ratios, «, so that 

„ioooo ^10, 

we may, since we find that approximately 

^3010 ^ 2, 

say that 0.3010 is the logarithm of two. Briggs was an ardent computer, 
and from 1614 to 1 617 he calculated to fourteen decimals the logarithms of the 
integers from 1 to 1000. In 1624 he published his "Arithmetica Logarith- 
mica" containing the logarithms to fourteen decimals of the integers from 1 to 
20,000, and of those from 90,000 to 101,000. This earliest of all logarithmic 
tables contains all the logarithms necessary for our modern five- and six-place 
tables. Briggs was still occupied in calculating the logarithms of the integers 
between 20,000 and 90,000, when in 1628 Adrian Vlacq of Holland published 
a table giving to ten decimals the logarithms of all integers from 1 to 100,000. 
These tables have been the basis of compilation of nearly all tables of loga- 
rithms published since that time. Naturally they contained some errors. 

It is noteworthy that these tables differ in arrangement from tables now 
in common use in that the logarithms are given to ten and fourteen decimals, 
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while the numbers, "arguments" of the table, contain only five figures. 
Modem four and five-place tables generally give the logarithm with only one 
more figure than is contained in the argument. That the latter arrangement is 
that best suited to practical numerical work appears from the following con- 
siderations : in such a table the difference of successive logarithms, that 
is, logarithms of arguments differing by one, is always less than forty- 
five units in the last decimal place, and for half of the values of the ar- 
gument, less than ten units, so that interpolation for a logarithm of a 
number containing the same number of figures as the logarithm in the 
table may be performed mentally.* Now if two numbers are to be 
multiplied by the use of logarithms, and if at least one of these numbers is not 
known certainly beyond m significant figures, the product cannot be found 
correctly to more than m significant figures,* and just that number of figures 
will be given by the use of logarithms to m decimals.! Hence for multiplica- 
tion of numbers of m figures it is desirable to use logarithms to m decimals, 
and consequently convenient, but not always necessary, to have a table of 
logarithms whose arguments contain m — 1 figures. 

One may then reasonably inquire whether tables like those of Briggs and 
Vlacq have, in the last figures of the logarithms, any value. To this may be 
answered first that these tables were doubtless intended to serve as a basis for 
the subsequent calculation of larger tables ; secondly that by interpolation the 
values of logarithms of intermediate numbers may be found from these tables 
to a large number of figures ; indeed, as I shall prove, t the value of the loga- 
rithm to m decimals of any number may be found from an m-place table 
(one in which logarithms are given to m decimals) by interpolation with first 
differences if the argument is given with ^m + 1 figures, t 

Briggs' method of calculating logarithms was extremely laborious, for at 
that time no developments in infinite series of the logarithm had been discov- 
ered. He extracted successive square roots of ten fifty-four times, obtaining 
the result, 

lOy = 1 + 12781 91493 20032 35 X lO"** = 1 + a, 

where y = ^^. 

♦ See Note A. 
t See Note B. 
^For a more exact statement, see Note B. 
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The right hand member of this equation may be written as one plus a decimal 
of thirty-two figures, fifteen ciphers preceeding the first significant figure of 
the fractional part. This i-esult amounts to the statement that the Briggsian 
logarithm of the second member of this equation is 

l/2« = 55511 15123 12578 2702 x lO"*^. 

Briggs discovered that the logarithms of numbers of the form 1 + x, where x 
is a decimal beginning with fifteen ciphers, are very nearly proportional to the 
decimal x. This ratio 

X 

approximately, for a small value of x, as we see from the series (1'). He ob- 
tained the value of the ratio by dividing l/2»* by ten raised to that power, and 
obtained a result to eighteen decimals of which the first sixteen agree with the 
correct value of M. His process was to extract successive square roots of the 
number whose logarithm was sought, until a root of the form 1 + * was ob- 
tained. The logarithm of this root was then found by multiplying x by the 
ratio 

log (1 + g) 
a ' 

and from the logarithm of the root was found at once the logarithm of the 
number sought.* 

Much interest was taken by mathematicians of the seventeenth and eight- 
eenth centuries in the calculation of logarithms. Most computers used methods 
of an arithmetical nature similar to that invented by Briggs. But as the 
calculus developed, less painful methods based upon the use of infinite series 
were invented and used by some scholars, among whom was Newton. There 
have been few tables made from new calculations since the publication of 
Vlacq's work, though many writers have computed the logarithms of some 
numbers. Indeed a new computation of logarithms already known cannot be 
regarded as a very useful service to the mathematical community. The most 
important calculations, after those of Briggs and Vlacq, were made by two 
Englishmen, Sang and Thomson, and by direction of the French government 
for the "Tables du cadastre." Sang published in 1871 a seven-place table of the 

♦ See Note C. 
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logarithms of the integers from 20,000 to 200,000, of which the last 100,000 
were freshly computed by him. Thomson calculated to twelve figures the 
logarithms of all iiitegers to 120,000. His work has been used to verify errors 
discovered in Vlacq's tables by a comparison with the "Tables du cadastre." 
An account of the latter tables may be of interest. A full description of them 
with explanations of the methods of their compilers is given by Lefort in an 
article in the fourth volume of the Annates de V Observatoire de Paris, 

In 1784 it was voted by the French authorities that new tables of the 
logarithms of numbers, of the trigonometric functions, and their logarithms 
should be calculated to correspond to the decimal division of the quadrant. 
The manuscripts of these tables, which have never been published, give to 
fourteen decimals the logarithms of the integers from 1 to 200,000, the nat- 
ural sines, and the logarithms of sines and tangents. The intention was to 
publish them as twelve-place tables, but even the twelfth figure is not reliable. 
The work was done under the direction of an engineer, Prony. His subordi- 
nates were divided into three sections : first, five or six mathematicians, in- 
cluding Legendre, who were occupied with the preparation of formulas and 
with other purely analytical work ; second, a group of seven or eight men 
who had at least enough mathematical knowledge to translate general formulas 
into numbers, and who did most of the calculating done from series ; finally, 
a group of seventy or more computers who were occupied chiefly with the 
work of interpolation. The work was performed wholly in duplicate by inde- 
pendent workers, and required over two years for its accomplishment. The 
two manuscripts are deposited in Paris, one at the Observatory, the other at 
the Institute. The tables received, as has been stated, the name of the "Tables 
du cadastre." The exact meaning of the word, "cadastre" is not clear. In 
modern French usage it means a public register of the ownership of real prop- 
erty. The term is applied to other registers, especially to those carried 
out in great detail, and it is presumably in this sense that it was applied to 
these mathematical tables. The chief interest of these tables lies in the facts 
that no other computation on so great a scale has ever been carried through, 
and that this is by far the greatest work of interpolation ever undertaken. 

The computation of the tables proceeded as follows. The logarithms of 
all integers from 1 to 10000 were computed to nineteen decimals, those of 
prime numbers from the series (2'), those of composite numbers by adding 
the logarithms of the factors. Thus was known the logarithm of every hun- 
dred from 1 to 1,000,000, and the logarithms of all integers between these 
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hundreds were obtained by interpolation. Only alternate logarithms were 
used as a basis for the interpolation, the others serving to check this work. 
The interpolation was not carried on by differences of known values, but the 
differences were themselves computed from series. The interpolation formulas 
used are of interest. 

Let Mo, Ml, Ma, etc., represent successive values of the function consid- 
ered. Thus, 

Mo = log n, 11^ = log(n + 1), etc. 

We represent the first differences by Amo, Am,, etc., so that 

Ml = Mo + ^^'H Mj = Ml + AMj. 

If the second differences are A^Mo, A^Mi etc., 

Amj = Amo + A^Mq, Am2 = Am, + Ahi^. 
In general, if A" be an with difference 

A"»-% = A^-iMp.i 4- A^Mp. 
By a combination of these formulas we find 

M, = u,+pAu,+Pi^ AX + ^^^ - V;^^ - '^ ASmo + ■ ■ ■ 

This formula is called b}' Lefort "Mouton's Interpolation Formula," but is an 
immediate consequence of Newton's formula.* 

In this formula the differences may be calculated from series, for 

Atto = log(n + 1) - log w = logfl + - j 

= ilf f- ^ 5— „ + r-r y 

\n Sn-^ 5*1* / 

A'^Mo = Amj — Amo 

= ^ (^TTi - 3(;rTip +■■■)- "^ G - i + ••• ) 

*Markoff, Differenzenreehnung, p. 16. 
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In this way may be obtained series for the differences of all orders. In the 
construction of the "Tables da cadastre" differences up to and including those 
of the sixth order were used. The tables, if completed, would have given the 
logarithms to twelve figures of all arguments with six figures, but as stated, 
the work was not carried beyond the finding of the logarithms of numbers up 
to 200,000. 

3. Practical Computation. Ifacomplete table of logarithms is to be 
calculated, there is doubtless no better method than that of interpolation, used 
in the Tables du cadastre. But if a single logarithm is to be calculated to a 
large number of decimals, the best method is one published by Weddle in 1845 
in The Mathematician, a method which is not directly based on infinite series. 
Before examining this method it is interesting to know why it should ever be 
necessary to compute a logarithm to a large number of decimals. This will 
be made clear by the statement of Gernerth in the preface to his excellent five- 
place tables.* Every logarithm is given by him correct to five decimals, 
and moreover it is indicated whether this value is larger or smaller than the 
correct value. His table was compiled from Vega's Thesaurus, published in 
Leipzig in 1794, which gives logarithms to ten decimals. But Vega's last 
decimal place is sometimes wrong, and Gernerth decided to consider this tenth 
figure as unreliable. Finding in the Thesaurus 

log 5.873 = 0.768 8600 008, 

Gernerth, fearing that the last figure might be nine units too large, was in 
doubt whether to write 

log 5.873 = 0.76886 + 
or 

log 5.873 = 0.76886 -, 

and was consequently obliged to compute the logarithm to ten places. It may 
be imagined that in some cases it would be necessary to compute a logarithm 
to even more than ten figures. 

Again, for a single calculation, the easiest method may sometimes be to 
compute the logarithms needed for the work and finally to compute the num- 
ber from the resulting logarithm. Thus, if it were required to find to thirty 
decimals the value of 1/23*", and only Gernerth's five-place tables were at our 

* A. Gernerth, Funfstilligt $reinet'n« Logarithmen, Vienna, 1901. 
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disposal, the best method would be to calculate to fifteen decimals log 23, and 
then from log 1/23'", to compute the number required. 

In most practical work it never becomes necessary to compute a logarithm. 
It is seldom necessary to use logarithms to more than" seven places, and gen- 
erally four or five places suffice. 

A description of Weddle's method, with an example of its application, is 
given by Gernerth, together with the small tables needed. This method I 
shall now describe and illustrate. 

We suppose, as we may without loss of generality, the number JV^, whose 
logarithm is sought, to be between one and ten, and write 

2f= a + d, 

where a is an integer less than ten, and d a decimal fraction. 

We now divide a + d by a + 1, writing the quotient, which is less than 
one, as a decimal carried to a number of places sufficient to insure the degree 
of accuracy desired." This is next multiplied by 1 + cv/lO*", where c, is the 
difference between nine and the first figure of the quotient not nine, which figure 
we suppose to be in the rth decimal place. Generally we should have r = 1. 
This product will have all figures, before the rth, nines and generally the rth 
also. This process of multiplication by factors of the form, 1 + 0^/10'', may 
be repeated until the final product is as near one as is desired. 

Then, with as small an error as is desired 



log iV^= log (a + 1) -2 log (l + j^,) • 



To follow this method it is necessary to have a table giving the loga- 
rithms of integers from one to ten, and the logarithms of numbers of the form, 
1 + Cr/lO*", c,. being always an integer less than ten. Such a table, giving 
the logarithms to fifteen places for all values of r up to and including sixteen, 
is given by Gernerth on a single page.f It remains to explain how this short 
table may be computed. 

* In the example, seventeen ; see Note D. 
tLoc. cit., p. 119. 
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Logarithms of numbers in the form 1 + -t^ may be easily found from the 

series (1'), when the value of M is known. The logarithms of the prime 
digits 2, 3, 5, and 7, and the value of M are ingeniously found by J. C. 
Adams in the following manner : 



Let 



a = 


10 
9 ~ 


■'( 


f 


lo)' 


b = 


24 ~ 


-'( 




looy' 


c = 


80 ~ 


'( 


f 
1 + 


8o)' 


d = 


50 

49 ~ " 


.i{ 


f 


looy' 


e = 


126 
125" 


t 




1000,/ 



The values of a, b, c, d, and e are easily obtained from the series (1). Then 
we shall have 

Z 2 = 7a - 26 + 3c, 

I 3 = 11a - 36 + 5c, 

? 5 = 16a - 46 + 7c, 

I 7 = J(39a - 106 + 17c - <Z) = 19a - 46 + 8c + e. 

The first of these formulas follows from the equation, easily verified, 

'10\V24\V81\* 



-(V")(ll)(l-J) 



The others may be similarly proved. From the two expressions for I 7 we 
obtain a check formula 

a — 26 + c = <Z 4- 2e. 
We have further 

1 



M 



= no = Z 2 + Z 5 = 23a - 66 + 10c. 
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Adams used these formulas to find to 260 decimals the natural logarithms of 
2, 3, 5, and 7, and the value of il/ to 282 decimals. It may be noticed that 
the discovery of such relations as those used by Adams would be a real appli- 
cation of simple indeterminate equations required to be solved in integers. 

We proceed to calculate by Weddle's method the logarithm of 5873, 
using the .-short tables mentioned given by Gernerth. Since Gernerth's tables 
are carried to fifteen places, we shall proceed with the multiplication until the 
first eight figures of the product are nines, when the remaining seven factors 
necessary for the computation may be found by inspection. The work of 
division and multiplication must be carried to seventeen figiires to insure the 
correctness of the fifteenth factor. To indicate multiplication by a factor 
1 -t- Cr/10'', we place this number to the right of the multiplicand, write 
underneath the significant figures of the product by cv/10'', correct to seven- 
teen decimals, and, adding this product to the multiplicand, obtain the product 
by 1 + c^/W. The factors obtained by inspection we shall denote by 
1 4- a,./10^ We have 

log 5873 = 3 + log 5.873 



5.873 
6 


0.97883 33333 33333 33 
- 1957 66666 66666 67 


2 

1 + 10^ 




0.99841 00000 00000 00 
99 84100 00000 00 


1 + 103 




0.99940 84100 00000 00 
49 97042 05000 00 


, 5 
1 + 10^ 




0.99990 81142 05000 00 
8 99917 30278 45 


9 

i + io» 




0.99999 81059 35278 45 
9999 98105 94 


i+io« 




0.99999 91059 33384 39 
7999 99284 75 


8 
1 + 10' 




0.99999 99059 32669 14 
899 99991 48 


9 

1 + lOa 



0.99999 99959 32660 62 

Then c^ = 0, c^ = 2, Cs= 1, c^ = 5, Cg = 9, Cg = 1, c, = 8, Cg = 9. The re- 
maining seven factors, 1 + a^/W, may be obtained by taking the differences 
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between nine and the successive iigures to the fifteenth. The fifteenth figure 
ia here correct because the errors of the appi*oximate work cannot be greater 
than five in the seventeenth place.* Then 

ag = 9, 09 = 4, oio = 0, au = 6, a,2 = 7, Ois — 3, a^^ = 3, ajg = 9. 

There are two methods of using Gernerth's short table, one taking the 
logarithms as they are given, the other adding or subtracting to each loga- 
rithm, according as it is too small or too large, one fourth of a unit in the 
fifteenth place. The plus and minus signs written after the logarithms and 
the additional plus sign are used in the second method. The logarithms of 
the thirteen factors different from one are written in order, each correct to 
fifteen decimals. 

0.00860 01717 61918 - (25-10-") 

0.00043 40774 79319 - 

0.00021 70929 72230 + 

0.00003 90847 44584 + 

04342 94265 - 

03474 35447 - 

00390 86502 - 

00017 37178 - 

26058 - 

03040 + 

00130 + 

00013 + 

4- 

+ 

0.00929 12495 40788 - (50-10-") 
Now, log 6 = 0.77815 12503 83644 - 

Subtracting, 

iog6-i:iog(i+^,)-|;iog(i+^) 

= log 5.873 = 0.76886 00008 42856. 
•Note D. 
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This value is liable to an error less than 7.55 units in the fifteenth deci- 
mal place if Gernerth's table be used by the first method ; to an error less 
than 3.85 units in the fifteenth decimal place if the second method is used. It 
is only a coincidence that the two methods give the same result. The second 
method enables us to say that this value is correct to fourteen decimals, but not 
that the value of this correct to fourteen decimals, 

0.76886 00008 4286, 

is the true value correct to fourteen decimals. Either method enables us to 
say that the true value correct to fourteen decimals is either that just given or 
diflers from that only in having the foui*teenth figure a five. 

4. The Limits of Error in Several Approximations. In pre- 
paring this sketch of the history of the development of logarithms and of the 
methods of their computation, I have frequently had occasion to question to 
what extent certain processes employed can be relied upon to give correct 
results. In attempting to answer the questions raised I have learned first, that 
we can in all cases assign to the error committed in the approximation a limit 
which cannot be surpassed, which I shall speak of as the limit of error ; second, 
that these limits are often smaller than I had before supposed ; and thirdly', in 
seeking to make the limit of error as small as possible, I found myself led 
in Weddle's method of computing a logarithm to very definite rules for com- 
putation. These results have interested me so much that I venture to hope 
that others may also find them interesting and not unprofitable. For a know- 
ledge of the limit of error in computation will not only give to the computer 
more confidence in his result, but will also, it seems to me, give him a much 
more truly scientific habit of mind. I have divided these studies into four 
notes referred to in the previous pages. A few preliminaries will be useful. 

It is customary in giving the approximate value x', of a number x, to 
give x' to a stated number of decimals m, so that |x — a;' | is as small as possi- 
ble. Then x' is said to be the value of x "correct to m decimals." Clearly 
we shall always have 

l^-^'I^ITIO^- 
If X and y are any two numbers such that 

l*-y|^2:io=' 
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we shall find it convenient to say that y ia equal to x correct to m decimals. 
But we must note that if as' and y' are the values of x and y correct to m deci- 
mals, we may have 

J_ 

lO™* 



i»'-yi- 



A few examples will make these points clearer. Let x — 0.43429 ; then 
if m is three, «' = 0.434 ; if m is four x' = 0.4343 ; if y = 0.4346, y is equal 
to 35 correct to three decimals, but for to equal three, y' — 0.435 and 

y - x' = 1/103. 

Finally if x and y are two numbers such that 



\ 
10 
then also 

\ 
10 



X — y\^ — — 






For example, taking again m as three, and x = 0.43429, if \x — y\ ^ 
1/10', y lies between, or is equal to one of the values 

0.43329, 0.43529, 

and y' is equal to one of the values 

0.433, 0.434, 0.435, 

no one of which differs from x' by more than 1/10*. 

We shall have occasion to use the developments for small values of x 

log(l + x) = M^x-j + ^^ 4 • • •}' 

< X < 1, 



10g(l-X)= Jt/Jx + y+g- +^- • -J. 



From these developments we may infer that, for positive values of x less 
than one, 

< log (1 + x) < Mx, 

< - log(l - x) < J—— . 



16 WHTTTEMORB [October 

Note A. 1. If two numbers, of which one is not known beyond m 
significant figures are multiplied, the product cannot be known beyond m 
significant figures. 

2. In an m-place logarithm table whose arguments contain m — 1 figures 
successive tabulated values differ by not more than forty-five units in the last 
decimal place, and for one half of the arguments by less than ten units. 

Proof of 1 . 

We may, without loss of generality, suppose both of the numbers to lie 
between 1 and 0.1, so that each number has its first significant figure in the 
first decimal place. If the two numbers are x and y, and they are known to 
in and n figures respectively, and if the approximately correct values be x' and 
y', then 

a;' = a; + e, y' = y j^ r), 
where 

l^l^l^lo^i' I''! = 2:10--' ** = '"• 

Then is 

x'y' =z xy + ey + TJX + erj, 

\x'y' -xy\ = \€y + vx + er)\. 

Since nothing is known of the signs of e and 17, we cannot know that the 
second member of the last equation is less than | «y | , which itself cannot be 
known to be less than or equal to l/C2-10'"+i), since y is greater than 0.1 ; 
hence it cannot be known that the value of xfy' is that of xy correct to more 
than m figures. It may also be shown that the upper limit of | x'y' — xy\ is 
(5.75)/10'", but that in more than two thirds of all possible cases this differ- 
ence is less than I/IC". 

Proof of 2. 

Let a and a + 1 be two successive values of the argument, and let the 
tabulated values of their logarithms, each correct to m decimals be log' a and 
log' (a + 1) . Then 

lloga-log'^l^^:^,, |log(a+l)-log'(a + l)|^2;^, a S lO-"-'. 

The tabular difference. A, is 

A = log'(a + 1) - log' a = log (a + 1) _ log o + ^, 1^1 ^ j^ • 
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Now is 

log(a+l)-loga = log(l+i)<f<3^. 

Then 

* 44 ,-, 45 
A<ioS+l^l<IOS- 

M 44 

If the first digit of a is five or greater, — is less than , , and 

, 8.8 , „, 

9 
< • 

lO- 

Note B. The Limit of Error in Interpolation by First Dif- 
ferences in a Table of Logarithms. Let us consider a table giving the 
logarithms correct to m decimals of all integers between 10" and 10"+^, a 
table then with arguments of n + 1 figures. Let a be any value of the argu- 
ment, and let log* o be the tabulated value of its logarithm. The formula for 
interpolation by first differences is 

(3) log (a + «) = log* a + » [log' (a + 1) - log'a] , 

where » is any positive number less tiian one. We proceed to investigate the 
limit of error of this formula. Let us write 

(4) log(a + a;) = log a- + (« + c) [log(a + 1) - log a] . 
We have from (1) 

log(a + «) — logo = ilfZ ( 1 + - j 

log(o + 1) - loga = J/ [1 - ^ + ^3 . . .], 

M 
whence, from (4), cancelling — , 



1^ 2a "^ 3a* ■ ■ J ~ 



»(1 — «) x(l — «*) 
2a 3a» 

x{\ — a;) 
2a ' 
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since each term of the series in the second member of the last equation is less 
in absolute value than the preceding.* Moreover the coeflScient of e in the 

first member is greater than 1 — ^ . Hence we find 

x(l — x) 

The maximum value of a;(l — x) is J. Then 

1 



< €< 



4(2a- 1) 

Since we know that 

M 
< log(a + 1) — log a < — , 

we may write 

log(a + x) =loga + X [log(a + 1 ) — log a] + e', 

where < e' < -r-r^ r.- • 

4a(2a — 1) 

If we write 

log a = log' a + €i, log(o + 1) = log'(a + 1) + eg 

^eknow |ei|^^-^^, |eg|^^-^^. 

Now 

log (a + x)= log' a + X [log'(rt + 1) - log' a] + €i(l - x) + Cj x + e'. 

But 

|ei(l — x)+ e-i x\ S |6i| (1 — x)+ |€2|x ^ 



2.10" 



If now e' is less than , log (a + x) will diflTer from the value ob- 

tained from (3) by less than 1/10"*. This will be the case if 

M 1 

< 



4a(2a- 1) 2.10'» ' 



* The ratio of the nth term to the preceding Is 



1 — X n 1 n= 1 

- < -2 T • - < 1, 



I x"~^ n + 1 a n^ — 1 a 

since x < 1 and a > 10. 
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or if 



or finally, if 



(4o- 1)2 > 4.10'". M+l, 
V4.10'».iJ/ + 1 + 1 



a > 



If m is even it may be easily seen that this inequality is satisfied by a S 10 2 ; 

ni — 1 

if m is odd, and equal to five or greater, it is satisfied by a S 1.05 x 10 z . 
Then in a table where m is even, if the number of figures in the argument, 
n + 1, is equal to m/2 + 1, the logarithm of any number may be found by 
interpolation with the formula (3) so that the error is not more than one in the 
last decimal place. Thus if m be four, there are needed for this degree of accu- 
racy arguments of three figures. If m is odd and as great as five the number 
of figures needed in the argument to ensure this degree of accuracy is 
(m + l)/2 if the first figures of a are as great as 105 ; for values of a beginning 
with figures^ less than 105, the arguments must be given with an additional 
figure.* 

There results from these considerations this remarkable fact : a logarithm 
obtained by the formula (3) from m-place tables with arguments of m/2 + 1 
or (»t + l)/2 figures,! when taken correct to m decimals, cannot differ by more 
than one unit in the mth decimal place from the true value of the logarithm 
correct to m decimals. This limit of error is the same that we should have if 
the process of inteipola.tion were exact. For in that case, since each tabulated 

* Gernertb, in his five-place tables, already referred to, uses arguments of foor figures, 
with an additional table for arguments of five figures to 10800. This is, according to the 
theory here presented, one figure more in all cases than is necessary for reliable-interpol&tton, 
but his arrangement evidently makes the interpolation mnch less laborious. 

The seven-place tables of Dletriclikeit (Berlin, 190S) have arguments of four figures. 

The eight-place tables published by direction of the French Ministry of War in 1891, and 
the eight-place tables of Mendiz&bal Tamborrel (Paris, 1891) both have arguments of five 
figures. According to this theory interpolation by first differences is reliable in these eight 
place tabie», and also in the seven-place tables vrhen a is as great as 1060. 

The ten-place tables of Vega have arguments of fire figures, so that Interpolation by 
first differences is not reliable. 

A brief treatment of the limit of error in a special case covered by my discussion is given 
by Markoff, Differenzenreehnung, pp. 36-39. 

A rough discussion leading to an apparently similar result is given on page 68 of the 
tables of Dietrichkeit. 

w— 1 

t According as m is even or odd, and when odd If a > 1.06 X 10 * ' 



20 WHITTEMORE [October 

value is liable to an error of 1/(2.10'"), the value obtained for \og{a + x) 
would be liable to an equal error, and consequently this value correct to m 
decimals might differ by one unit in the with place from the true value correct 
to m decunals.* 

We have, finally, supposing our m-place table to have arguments of 
m/2 + 1 or (m + l)/2 figures.t to see with whatlimitof error a number iVT, may 
be found from its logarithm. We may, without restriction, suppose the 
characteristic of log iVto be n, so that JVlies between two successive argu- 
ments of the table, a and a + 1. Then a positive value of x, less than one, 
may be found from the equation 

log iV"= log'a + X [log'(a + 1) — log* a]. 

The determination for iV is then 

If' = a + X. 
Now the value determined by (3) for log JV' is precisely log JT. Hence 

1 M 

|logi^-logi\^V2:iO^+ 4^(2a-l) - 

If we write 2^= A-10"+^ and JST' = A'-lO''+\ A and A' both lie between 1 
and 0.1, and 

iV 

Let us suppose that A' is less than A. Should this not be the case the iden- 
tical result would be obtained by interchanging these quantities in the follow- 
ing lines. We have 

l^^--^, e > 0. 



^ \iogjff-\ogif'\ = \ii^-m'\ =1 



A' 



Then is 



A' = Ae-' = a(i - e + ^ 



* If in (3), a + X is tlie valae of a number N, correct to m figures, the limit of error in 
taking the second member of (3), correct to m places, for the value of log N, correct to m 
places, is three in the last place even if a contains m — 1 figures. But if o + x is the value of 
JV, correct to )?i + 1 figures, the limit of error is, as before, one in the mth place, if n = m/2 
when m is even, or if n = (m — 1 )/2 and a > 1.5 X 10" when m is odd and at least equal to five. 
Then in logarithmic work where numbers are only approximately known it is necessary, to ob- 
tain the best results, to use these numbers correct to one more figure than is contained in the 
tabulated logarithms. 

i n — 1 

t According as m is even or odd, and when odd if a > 1.05 X 10 i ' 
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and 

1-4' — ^1 < -^e < A\ --■-■, -^ + -r—r^ 7^1- 

' ' (2Jf.lO'» 4a(2a— 1)J 

Since ^^ < 1.152 and 4a(2a — 1) > 7a^, we have 

Now, however n is chosen, we cannot infer that | ^' — -4 1 S I/IO", so that 
it will be impossible to assert that the values of A and A', each correct to m 
decimals, diflfer by less than two in the last place. But if n is chosen equal to 
{m — l)/2 or n»/2, according as m is odd or even, it is easy to see by ex- 
amining the different possibilities, that | -4' — -4 1 < 2/10'", and conse- 
quently that the values of A and A\ each correct to m decimals, cannot 
differ by more than two in the last place. The truth of this statement is im- 
mediately evident when m is even. If m is odd suppose, for example, that 
the first digit of a is 3 ; then is 

, ., ,, - ,f 1-152 10 ] .53 

1^ — ^1 < "•*lT05r + g3:i05'J ^ lo^s • 

In fact it appears in this way that, if the first digit of a is seven or less, 
1^' — ^1 < 1/10"». We may then say, returning to the numbers iVand^', 
that the values of these two numbers each correct to m significant figures can- 
not differ by more than two in the last figure in any case, nor by more than 
one in the last figure if the first figure is less than eight. 

Note C. Tlie Limit of Error in Briggs' Method of Calculating 
a Logarithm. Briggs' method is based on the repeated extraction of square 
roots to a large number of decimals. It will perhaps be of interest to prove 
the fact, notalTrays realized, that the square root of a number larger than one 
may be always found correct to at least as many decimals as the number is 
given, sometimes to more. 

Let x' be the given value, correct to m decimals, of a number x. Then 

i« t*-^'l^2305.- 

We seek the square root of x by finding the square root of x'. If this pro- 
cess is carried out correct to n decimals, and we write the resulting number 
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as ^x", we have 

Now, the absolute value of the error, 

I V^ - VP| S Iv^ - V'x'l +1 V^ - V^'l 

^ ^ 1 1 

- 2.10«» ■ v'»+ V^"^ 2.10"' 
since 

, 1 1 



\\Jx — \ja^\ = 



2.10"' yjx + V^' 



Then if x is greater than one, we have 

\[x + yfx' > 2, and if n = wi + 1, 

But the values of y'a; and y/a;", each correct to m decimals, may differ by one in 
the »ith decimal place. We may notice that if x > 10*p, we may, if we choose 
n = TO + J) + 1, assert 

that \\lx-\l¥'\< 2.\Q^+p ' 

so that the values of ^x and \jx", each correct to to + ^ decimals, differ by 
not more than one in the last decimal place. 

It has been explained that Briggs' method of computing the logarithm of 
a number, N, consisted in finding a root 

V^ = 1 + «'. 

where p is of the form 2*, and x a decimal beginning with at least fifteen ciphers, 
then in writing 

log(l + a;) = ^log(l + «), 

where 1 + a = 10", where y = 1/2** 

so that a = 12781 91493 20032 35 X 10-^. 

Let us write x = na; since x is to begin with fifteen ciphers, n must be less 
than eight. Briggs' formula becomes 
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log(l + na) =. n log(l + a) = log(l + a)». 
By Taylor's Theorem, we have 

(H-a)» = l + na+ "(** - ^)^' (i + <?a)»-2, 1 > <> > 0, 

whence 

1 + na = (l + o)"(l — e), 

where (1 + a)"e = **(" - ^)«' (i ^. ^o)—*. 

Then log(l + no) = n log(l + a) + log(l — e). 

Mt 
Now I log(l — e) I < j^^-;* 

and since n is less than eight, and ^—^^ ^—j — < 1, we have e < 28 cr*. 

Substituting the value of a, we find e < 46 X 10"**, whence 

|log(l-e;|<2.10-» 

Then for all values of n less then eight, 

I log(l ■if.na)-n log(l + a) | < 2.10-n. 

Now we may suppose that ^is always between one and ten ; then to find a 
root of the form 1 + «, it will never be necessary to take h greater than 54. 

Now is log i\r=j> log(l 4- x) = 2* log(l + «), 

and the error in log N, computed by Briggs' method, will not be gi-eater than 

2« X 2.10-" < ^ 



2.10" 

Hence the inaccuracy due to the method will not prevent the result from being 
correct to fourteen decimals. 

Note D. The Limit of Exror in Weddle's Method of Caloulat- 
ing a Logarithm. We suppoie that the number N, whose logarithm is to 
be calculated, lies between one and ten. We write ^= a + rf, where a is an 
integer less than ten, and d a decimal fraction. We divide a + <2 by a + 1, 
unless a -^ d begins with the figures 1.0, when the method may be advan- 
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tageously modified by dividing by 1.1 instead of by two. The quotient N' 
is then multiplied by I factors of the form 1 + Cr/lO*", so that the product has 
its first k decimal figures equal to nines. Representing the product of I factors 
of the form 1 + Cr/lC by Pi, we have 

7\7- " _ A7-> p _ 1 ^*+l ""k + i 

~ ' - To*+T " 10*+* 
It is always possible to obtain this result where I ^k + Z, and if Z = ^ + 3, 
we shall have a^^i = 0. The truth of this statement appears as follows : We 
have always iV' > . 55 ; in the least favorable case multiplication by four factors 
of the form 1 + cv/10'' will raise to nines the first two figures. When the 
first two or more figures in the product are nines it may be necessary to 
multiply b}' two factors to raise to nine the first figure of the product different 
from nine, but in the least favorable case this process will raise to nine this 
figure and the one next following. Generally, each multiplication raises to 
nine one figure of the product, so that I = k, but we may assert that in all cases, 
Z 5 i • + 3. An example will make these explanations clearer. Suppose 

N' =0.735 ci = 2 

IT'P^ = 0.8820 ci = 1 

iV"'P2 = 0.9702 C2 = 2 

iVr'Pg = 0.98960 4 4 = 1 

JST'Pi = 0.99950 004 Cg = C4 = 4 

JV^'Pb = 0.99989 98400 16 c^ = 1 
iV'Pg = 0.99999 98.S00 00 

In this rather unfavorable case if A; is two, we have to take Z = 4 = A; + 2 ; 
to make the fourth decimal a nine two factors, 1 + c^/10*, are used but with 
the fourth not only the fifth but the sixth figures become equal to nine, so 
that if k is six, we have I = 6 — k. 
Having now 



10*+i io*+2 



we write 



\ ^ io*+V V ^ io*+v 

and for x we find 



= 1 



+ 



jL^ 10'+> ^ 10<+.'+'» 
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where k+l^i^J^ni''', 

and where not more than two indices are equal in any term of the sum. Since 
DO value of a is greater than nine, we have 

SfljO, 81 •^ 1 _ 1 

where a and /3 take independently all positive integral values. Similarly, 
since m is at least equal to A; + 2, 

SaiOjO^ 729 -^ 1 _ 1 

So that, finally, we have 

1 f 1 1_ I 

102*t "•■ iP+i "^ l02»+2" ■ ■) 

1 1 1.01 .- , , 

< IQiS i — ^lO**' u fc> I. 



X < 



10*+i 



Now log If" 4- J ^°S (l+^)= log(l - *) • 

rait + l 

But |log(l_«)|<^__<_^, 

since x is known to be less than 1.01 X 10-**. 

Then logiVT" = - J) log (l + ^,)+ e, |e[< ^ 

and logiyr' = logiyr'+ 2 log (i + j§-,) ; 

and since log Jf= log (a + 1) — log .y, we have 

logi<7=log(a + l)-2log(l + j|,)- 2) log(l+^r) + ^- 



r«s*+l 



If in the last sum we neglect all terms for which r is greater than m, 
we commit the error 

r!«i» + l ' rvflt + l 
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We compute log iVfrom the formula 

(4) log N = \og'{a + 1)- 2 I«g' (l + w) - 2 '«g' (l + w) ' 

where log' means, as before, the tabulated value of a logarithm. If the tables 
used are correct to p places, each logarithm used is liable to an error of 

1/(2-10''), and we use 

\ j^ I j^ m — h, 

tabulated values, or less if some values of a^ are zero. Now I is not greater 
than A; + 3, hence not more than m + 4 logarithms are used. Then the limit 
of error in (4) is not greater than 

F-MJ_ . 1 I ^ + ^' 

2 1 10«* "^ 10™ "^ lOP 

If the number of logarithms actually used is less than m + 4, this number 
may replace m -f 4 in J?. 

"We shall next consider what values may best be assigned to I- and m, 
when p is given. It is desirable to make E as small as possible, and, in order 
to shorten the labor of computation, to take h no larger than is necessar}'. 
The least value possible for 

1 m 

10^ + 10^ 

is obtained by writing m = p. Then 

1 m _p + 1 

10" "^ 10^ "" ~1^' 

For clearly a value of m Sjp + 1 would give the expression a greater value. 
Whereas, if m = p — h. 

1 m _ 10* + j9 — A p + I 

lO^ '^ Wp 10? ^ IQP 

We shall then do best to put m — p, and 

^_M 1 I i^ + 51 
2 1 10** ^ \0P ]' 

To choose 2k less than p would make E unnecessarily large, and since E has 
inevitably a considerable value in the ^th decimal place there is no gain in 
taking 2A; greater than ^ + 1. Then if j? is even, let k = pl2, and 

li> + 6 



j; = 



2 \QP 
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if p is odd let A; = (j? + l)/2, then 

1 p+ 5.1 

2 lOi' 

In all cases if the number of logarithms actually used in (4) is less than J9 + 4, 

the difference may be subtracted from the numei-ator of E. In the example 

worked out in §3, we had ^ = 15, hence took /» = 15 and h — %. There 

were used in (4) fourteen logarithms. Hence we subtract five from the numer- 

1 15.1 
ator of E and E = - -j^ and the error in the fifteenth decimal of log 5873, 

computed by the first method is less than 7.55. The value found, if taken 
correct to fourteen decimals cannot differ by more than one in the fourteenth 
place from the true value correct to fourteen decimals. 

Xow the chief part of E comes from possible inaccuracy in the tabulated 
logarithms. This part may be cut in half by a slight modification suggested 
and used by Gernerth. With each logarithm tabulated it is indicated whether 
the given value is smaller or larger than the true value. If the tabulated 
value, V, is too small, the true value, ?, is between V and V + 1/(2.10">), 
and if in place of V we use V + 1/(4- 10^'), the difference from I cannot be 
more than 1/(4-10^'). Similarly, if the tabulated value is too large, and in 
place of V we use V — 1/(4-10^^), the difference from I is not more than 
1/(4-10''). If this modification is practised, we have 



„_lf_2_ _2_ m + 4 ) 
41 lO^*"^ lO" ^ 10^ J' 



E' being now the limit of error in (4). 

It appears that m should be given the value p -\- \. The last logarithm 

in the sum, log'f 1 + 17^^ ) . ^i'l be tabulated as zero, for 

^°K' + ^)< ^«K' +irO -" il W.' 

but in the modified calculation, its value must be written as 1/(4-10''). As 

before, Tc should be chosen as (^ + l)/2 orp/2, according as^ is odd or even. 

If p is even, 

l p + 1.2 ^ 

4 10^ ' 
ifpisodd, E'= JP + ^•^. 

-'^ ' 4 IQP 
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If the number of logarithms actually used is less than p + 5, the difference 
may be subtracted from the numerator of U'. In the example of §3, using this 

second method, fifteen logarithms are used, so that JS' = -77555 • But even this 

result does not enable us to say that the value of log 5873, computed by the 
second method, and taken correct to fourteen decimals is equal to the true 
value correct to fourteen decimals. 

It remains to consider one more point. In order to determine the values 
of tti, ffli+i, • • • dp, to how many places must the division of JPiThy a + 1, 
and the multiplication of the quotient by the factors 1 + c^/lO'' be carried? 

Suppose all this work to be carried out correctly to n decimals. Then the 
error in the value used for iV' is not greater that 1/(2.10"). This error is 
multiplied by c^/lO'' and the product taken correct to n decimals, so that, in 

the value used for JV'(cv/10'"), the limit of error is ^ -iipA ^ + fT^j* Then in 
the value found to n places for iV' ( 1 + tk^J the limit of en'or is 



is 



:^.^^0 + ^)} 



2. 
Repetition of the reasoning shows us that the limit of error in jy'P, is 



rio^{i+(i+w.>-" + ^'} 



2.10" 



where 1 + r-^, is the last factor in Pj. We have alway Pj less than two, since 
iT' > .65 and iV'P; < 1 ; generally Pi is very much less than two, and since 
I <k + S < -^-i- + 3, the eiTor in JSf'Pj will be less than 5^,^^. It will 
then suffice generally to choose n = ^ + 2 to have a correct value of Op. In the 
example of §3, I is 7, and Pi < — , so that the error is less than 

8 ^^ 

^•y 4.5 



< 



2.10" 10" 

To have the fifteenth place correct it is here sufficient to take w = 17. 
Hasvars Univehbity, 

CAMBBIDaX, MaBS., AFKO., 1907. 



